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Abstract 

Recent experiments on two-frequency parametrically excited surface waves produce an 
intriguing "superlattice" wave pattern near a codimension-two bifurcation point where both 
subharmonic and harmonic waves onset simultaneously, but with different spatial wavenumbers. 
The superlattice pattern is synchronous with the forcing, spatially periodic on a large hexagonal 
lattice, and exhibits small-scale triangular structure. Similar patterns have been shown to 
exist as primary solution branches of a generic 12-dimensional Dg-i-T^-equivariant bifurcation 
problem, and may be stable if the nonlinear coefficients of the bifurcation problem satisfy 
certain inequalities Q. Ifere we use the spatial and temporal symmetries of the problem to 
argue that weakly damped harmonic waves may be critical to understanding the stabilization of 
this pattern in the Faraday system. We illustrate this mechanism by considering the equations 
developed by Zhang and Vinals for small amplitude, weakly damped surface waves on a 
semi-infinite fluid layer. We compute the relevant nonlinear coefBcients in the bifurcation 
equations describing the onset of patterns for excitation frequency ratios of 2/3 and 6/7. For 
the 2/3 case, we show that there is a fundamental difference in the pattern selection problems 
for subharmonic and harmonic instabilities near the codimension-two point. Also, we find that 
the 6/7 case is significantly different from the 2/3 case due to the presence of additional weakly 
damped harmonic modes. These additional harmonic modes can result in a stabilization of the 
superpatterns. 



*To appear in a special issue of Physica D dedicated to the memory of John David Crawford. 



1 Introduction 



Faraday waves are parametrically excited on the free surface of a fluid layer when it is subjected to 
a vertical vibration of sufficient strength. This pattern-forminff hydrodynamic system has proven 
to be an especially versatile one in laboratory experiments ^, exhibiting the common patterns 
familiar from convection (stripes, squares, hexagons, spirals), as well as more exotic patterns such 
as triangles Q, quasipatterns [0, (zj, ^, s uperlattice patterns |jl|, ^, time-dependent rhombic 
patterns jllj and localized waves |l9|, |f2|]. See for a recent review paper on Faraday wave 
pattern formation. 

The temporal period of the Faraday waves is typically twice that of the vibration in the case of 
purely sinusoidal forcing. The observation of this subharmonic response is attributed to Faraday 
and was first explained theoretically by Benjamin and Ursell's linear stability analysis for inviscid, 
potential flow |p^ . More recently it has been shown that waves, synchronous with the forcing, 
can be excited in thin layers of fluid vibrated at low frequency jl6[ |l^ |l8|; in certain viscoelastic 
fluids 1^ ; and in fluids forced periodically, but with more than one frequency component |Q, . 
In each of these Faraday systems it is possible to tune the forcing parameters in order to access 
the transition between subharmonic and harmonic response. At this codimension-two point, both 
instabilities set in simultaneously, but with different spatial wavenumbers. 

Many of the experimental [|, 0, |ll|, ^ and theoretical studies ||, |2|, |24[ H of 



exotic patterns in the Faraday system attribute their formation near the codimension-two (or 
"bicritical" ) point to resonant triad interactions involving the critical or near-critical modes with 
different spatial wavenumbers. In particular, the focus has been on spatial triads ki, k2 and 
ka = ki±k2, where |ki| = |k2| is the wavenumber of one critical mode, and jkaj is the wavenumber 
of the other critical mode. The angle Or, which separates ki and k2, is readily tuned by changing 
the frequency components mu and nuj of a two-frequency periodic forcing function. It has been 
suggested, for example, that by tuning this angle, different types of exotic wave patterns may be 
selected . Such a simple mechanism for nonlinear pattern selection, which is based on examining 
the linear instabilities of the spatially homogeneous state, is naturally attractive, but warrants 
careful examination as we show. 

Silber and Skeldon ||2^ recently showed that whether or not resonant triads associated with the 
bicritical point affect pattern selection depends on the temporal characteristics of the competing 
instabilities. For instance, the bicritical point of laboratory experiments typically involves a sub- 
harmonic mode (Floquct multiplier —1) and a harmonic mode (Floquct multiplier +1). On the 
subharmonic side of the bicritical point, the onset pattern selection problem is strongly influenced 
by the presence of the weakly damped harmonic mode. In contrast, on the harmonic side, the 
onset pattern selection problem is completely insensitive to the presence of near critical subhar- 
monic modes. These general ideas were demonstrated in [ p6[ through a bifurcation analysis of a 
hydrodynamic model of one-dimensional Faraday waves. 

Here, we extend the bifurcation analysis in p6| to two-dimensional spatially-periodic pat- 
terns and to higher forcing frequencies within the two-frequency forcing function. With the 
experimentally-relevant higher forcing frequencies {e.g. 6uj and 7uj) employed in this paper, we 
find the new possibility that spatially-resonant triads involving nearly critical harmonic modes 
may influence the harmonic wave pattern selection problem. This is not an option for the lower 
forcing frequencies (e.g. 1uj/2uj and 2a;/3aj) used in previous weakly nonlinear analyses of the 
two-frequency Faraday problem p3, E6[ . 
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We follow J.D. Crawford's seminal work on Faraday waves ||2^, ^ |2|, ^ by posing the pat- 
tern selection problem in terms of a symmetry-breaking bifurcation of the trivial fixed-point of a 
stroboscopic map. By restricting solutions to those that are spatially-periodic on some hexagonal 
lattice we obtain a finite-dimensional bifurcation problem that can be analyzed using the meth- 
ods of equivariant bifurcation theory pl| . For a review of this approach to hydrodynamic pattern 
formation problems, see Crawford and Knobloch p2) . 

This formulation of the bifurcation problem allows us to address recent two-frequency Fara- 
day wave experimental observations of a transition between simple hexagons and the triangular 
superlattice wave pattern depicted in Figure |l|a. Specifically, we follow ^ and consider a bifur- 
cation problem that is equivariant with respect to a twelve-dimensional irreducible representation 
of De-i-T^, which is analyzed in The observed harmonic wave states correspond to pri- 

mary transcritical branches of the generic bifurcation problem. In order for the observed hexagon- 
superlattice pattern transition to be reproduced by the bifurcation problem, we must consider a 
degenerate case in which the quadratic coefficient vanishes. Moreover, the cubic coefficients must 
satisfy certain inequalities, e.g. certain combinations of nonlinear cross-coupling coefficients must 
be small compared to the cubic self-coupling coefficient. 

In this paper we compute the quadratic and cubic nonlinear coefficients in the bifurcation 
problem from the Zhang- Vinals equations which apply to deep layers of low viscosity fluids 
subjected to a periodic acceleration. We show that the necessary inequalities for stable superlattice 
patterns can be satisfied for the forcing frequencies employed in the experiments {Qlo/Ilo), and 
that a resonant triad involving a weakly damped harmonic mode plays a key role in stabilizing 
the superpattern. Specifically, we find that the presence of a near critical harmonic mode leads 
to a cancellation in one of the cubic cross-coupling coefficients, causing this coefficient to become 
small in magnitude as required. This selects a preferred angle Or for the superlattice patterns. In 
other words, it suggests which of the countably infinite 12-dimensional irreducible representations 
of Dg-i-T^ is most pertinent to this Faraday wave problem. 



The paper is organized as follows. Section 2.1 presents background linear stability results for 



the two-frequency Faraday experiment, while section p.2| reviews results from p6| on the influence 
of spatio-temporally resonant triads on pattern selection. Section 2.3 then formulates the generic 
bifurcation problem relevant to our investigation. The bifurcation results derived from the two- 
frequency Faraday problem modeled by the Zhang- Vihals equations are presented in Section ^; 
the coefficients of the leading nonlinear terms are evaluated numerically from expressions derived 
perturbatively in the Appendix. We consider two different cases. In Section 3.2 we consider an 
example involving forcing frequencies in ratio m/n — 2/3, focusing on differences between the 
pattern selection problems for subharmonic and harmonic wave onset in a vicinity of the bicritical 
point. Section IsTs] then turns to an example involving higher forcing frequencies in ratio m/n ~ 6/7, 
and shows how weakly damped harmonic modes can stabilize harmonic wave superpatterns involving 
the angle Or associated with a harmonic wave resonant triad. Finally, Section ^ concludes the paper 
with a brief summary of our results and some discussion of issues we hope to address in the future. 



2 Background 
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2.1 Linear Results 



In the two- frequency Faraday wave problem a container of fluid is accelerated in the vertical direction 
with an excitation of the form 

5(*) = .90 + gz (cos(x) coa{niujt) + sin(x) cos{nujt + cj))) . (1) 

Here m and n are co-prime integers, so the forcing function is periodic with period T — and go 
is the usual gravitational acceleration. For small amplitude acceleration gz the surface of the fluid 
remains flat and the fluid layer is merely translated up and down with the drive. For higher values 
of gz waves are parametrically excited on the surface of the fluid layer. 

Besson, Edwards and Tuckerman starting with the Navier-Stokes equations for the free 
boundary problem, determined the linear stability of the flat surface in the case that the fluid layer 
has flnite depth but is unbounded horizontally. They used a Floquet-Fourier ansatz and solved 
the linear stability problem numerically to determine, for each spatial wavenumber k, the value of 
gz where a Floquct multiplier first crosses the unit circle. The resulting neutral stability curves 
show that the primary instability is to either subharmonic or harmonic waves depending on the 
value of X a-nd the values of m and n. (Harmonic/subharmonic response is relative to the forcing 
period T = 27r/w.) Typically, if x is small so that cos(x) cos(mcjt) is of greater significance than 
sm{x) cos{nujt + (jj), then the response is harmonic if m is even and subharmonic if m is odd. 
Similarly, if x is close to tt/2, the primary instability is (sub)harmonic if n is even (odd). At the so- 
called bicritical point, x — Xci both harmonic and subharmonic instabilities onset at the same value 
of the excitation amplitude, but with different wavenumbers. The harmonic superlattice pattern 
of Figure p^, observed by KudroUi, Pier and Gollub [|], was obtained near the bicritical point for 
m/n = 6/7 forcing in ([l]). The pertinent neutral stability curve, computed using the experimental 
fluid parameters, is given in Figure Clb. 



2.2 Spatio— Temporally Resonant Triads 

When the hydrodynamic problem is posed on a horizontally unbounded domain there is no pre- 
ferred direction (in the horizontal) so that each critical wavenumber from linear analysis actually 
corresponds to a circle of critical wavevectors. There are two such critical circles at the bicritical 
point, as shown in Figure |^. In this situation it has been argued that resonant triads may play a 
central role in the Faraday wave pattern selection problem |7[ ^, ^ |4|. Resonant triads are 
comprised of three critical wavevectors that sum to zero; two examples are shown in Figure ^. In 
the first example, -l-k^a — k„ = 0, and in the second example k„j — k„2 — k^ = 0. Here the to, n 
subscripts indicate that the critical wavenumbers can be roughly associated with the muj and ntu 
excitation terms in (|l|). We identify with each resonant triad an angle 6r G (0, which separates 
the critical wavevectors with the same length. For instance, the angle in Figure 0b satisfies 



while the angle in Figure Gb satisfies 
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(a) (b) 

Figure 1: (a) Blow up of the experimental superlattice Faraday wave pattern described in |l| 
(courtesy of Kudrolli, Pier and Gollub). The forcing function (||) has m/n = 6/7, x = 61° and 
(j) = 20°. Note that the pattern is periodic on a (large) hexagonal lattice, and that in each hexagonal 
'tile' there is small triangular structure, (b) The corresponding neutral stability curve, calculated 
from the full (linearized) hydrodynamic equations, for the experimental parameters reported in 
(Sub)harmonic resonance tongues are given by solid (dashed) lines. The neutral curves are com- 
puted using the method described in [pO[. 




Figure 2: (a) A plot of a neutral stability curve 5z(fc) showing minima at fc = km and k — kn- (b) 
An associated spatially resonant triad , and k„ = k^^ + k^a ■ (c) An associated spatially 
resonant triad k„j , k„2 and k^ = k„^ — k„2 • 
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To illustrate the potential for resonant triads to influence pattern formation in parametrically 
excited systems we consider a bifurcation problem involving the three critical Fourier modes associ- 
ated with the resonant triads of Figure ^. Much of this discussion is a review of the key theoretical 
ideas in ||2^ . Because of the periodic forcing of the system, it is natural to formulate the bifurca- 
tion problem in terms of a stroboscopic map Specifically, we denote the free surface height 
z = ft.(x, t) (x e R2) at time t = pT {p £ Z) by 

h{x,pT) ^ A{p)e'^'^-'' + B{p)e'^'2-'' + C{p)e''-^'^+^'^'>-'' + c.c. + ■ ■ ■ . (4) 

Here A, B and C are the complex amplitudes of the linear modes that are neutrally stable at the 
bicritical point and which form a resonant triad. In this discussion we assume that the angle 6^ 
between k/j^ and k/^ is not 7r/3 so that the critical modes interact nonlinearly to generate other 
modes on a rhombic (rather than hexagonal) lattice. These additional modes, denoted by • • • above, 
are linearly damped at the bicritical point. We may then use the spatial reflection and translation 
symmetries to determine the general form of the bifurcation equations that govern the dynamics 
on a center manifold. Specifically, to cubic order, the codimension-two bifurcation problem takes 
the form 

A -> (tA + oBC +{a\A\'^ + b\B\'^ +c\C\^)A 

B aB + aAC +{a\B\'^ + b\A\'^ +c\C\'^)B (5) 

C ^ /iC + MB + + e|Cnc , 

where A is the complex conjugate of A, and the coefficients are all real. The Floquet multipliers a 
and p are either -1-1 or —1 depending on whether the linear modes A, B, and C are harmonically 
or subharmonically excited, respectively. 

In deriving (|^) we considered only the spatial symmetries associated with the resonant triad. 
Following p^ , we enforce the temporal symmetry associated with the triad through a normal 
form transformation of (^). Specifically, there exists a near-identity nonlinear transformation that 
removes all nonlinear terms in (^) which do not commute with Lp" , where L is the Jacobian matrix 
associated with the linearized problem (see, for example, Crawford's review paper on bifurcation 
theory (35|). Here 

/ct 0\ 

i = (T , (6) 
VO p) 

where |(t| = = 1. The normal form symmetry may be interpreted in terms of time-translation. 
Specifically, advancing by one period in time maps period-doubled modes to their negatives, e.g. 
if /i = — 1, then advancing one period takes C —f —C. 

In the case that p, = +1 (cr = ±1), the bifurcation problem (|^) is already in normal form. 
This observation is trivial ifo' = +l. Ifcr^— 1, then the normal form symmetry is equivalent 
in action to that associated with the spatial translation symmetry x — > x + d, where d satisfies 

k;^ • d k;^ • d = IT. 

In contrast, in the case that p — —1, a normal form transformation removes the quadratic terms 
in the bifurcation problem (^. The normal form of the bifurcation problem, through cubic order, 
is then 

A aA+{a\A\^ + b\B\^ + c\C\^)A 
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B ^ t7B + (a|S|2 + 6|^|2 + c|C|2)B (7) 
C ^ -C+(d|A|2+d|B|2 + e|Cnc . 

We note that C = is a dynamically-invariant subspace of (^). This is true to all orders of the 
normal form since C = is the fixed point subspace of a (spatio-)temporal symmetry. Specifically, 
if (T = +1 then C = is the fixed point subspace associated with the time translation by one-period, 
i.e., {A,B,C) — > {A, B,—C). And if ct = —1, then C = is the fixed-point subspace associated 
with the spatio-temporal symmetry involving time translation by one period followed by spatial 
translation by d, where again k/j^ • d = k/^ • d = tt. 

We now examine ^ more closely in the case that /i = +1 so that we cannot remove the 
quadratic nonlinearities by normal form transformation. We focus on a detuning from the bicritical 
point such that the C mode is weakly damped, while the A, B modes are neutrally stable. In this 
case, \a\ — 1, /i < 1, we can further reduce the bifurcation problem to one involving the critical 
modes A and B, with C constrained to the center manifold: C — ^ . AB + • • •. We then obtain 
the reduced bifurcation problem 

A aA + a\A\^A + l3{er)\B\^A 

B aB + a\B\^B + l3{er)\A\'^B , (8) 

where the cross-coupling coefficient is 

(9) 

We see that in this case, the near critical spatio-temporally resonant mode C in (||) can contribute 
significantly to the cross-coupling coefficient P{6r) since < 1 - M < 1 in (I). For example, for 
fi sufficiently close to 1, the second term in (||) dominates and f3{9r) becomes large in magnitude. 
However, we also point out that if b and aS have opposite signs, then P{Or) could actually vanish 



for some /i — 1 > 0. Examples of these two very different situations are given in sections 3.2 and 3.3, 
respectively. 

We contrast the above with what happens when = — 1 at the bicritical point. In this case 
a = (5 = in the normal form and C = is an invariant subspace with associated dynamics of the 
form (H) with I3{9r) = b. In this case, the triad is spatially resonant, but not temporally resonant, 
and the cross-coupling coefficient is insensitive to any parameter proximity to the bicritical point. 

These observations about (3{9r) are important for understanding which patterns might be ob- 
servable near onset since branching direction and stability of patterns are determined by various 
nonlinear (cross-coupling) coefficients in the amplitude equations. We discuss this further at the 



end of Section 2.2 



Finally we note that similar results to the fi = —1 case above apply when there are weakly 
damped modes with complex Floquet multipliers. Specifically, these modes do not contribute sig- 
nificantly to the cubic cross-coupling coefficient /3{9), even when they are spatially resonant with the 
critical modes. Only damped modes with Floquet multiplier /i sufficiently close to -t-1 contribute. 
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2.3 Hexagonal Lattice Bifurcation Problem 

The analysis of the previous section led to certain conclusions about the nonlinear coefficients in 
the general rhombic lattice bifurcation problem 

Vl avi + {a\vif + P{e)\v2f)vi 

V2 <JV2 + ia\v2\^+P{9)\vi\^)v2. (10) 

Here ui, V2 are the complex amplitudes of two critical Fourier modes with wavevectors ki, k2 
(|ki| = |k2| = fee) that are separated by an angle 9 G (0, f] {9 7^ f-)- In particular, it follows 
from (^) that if a weakly damped harmonic mode is removed via center manifold reduction, then 
f3{9) may become large in magnitude when the spatial resonance condition is met, i.e. when 9 = 9r. 
This is in contrast to the situation where there are weakly damped subharmonic modes, which have 
no special influence on the pattern selection problem at onset. 

We now lay the framework for examining possible implications of these results for stability of har- 
monic hexagonal and triangular superpatterns. We follow Q and introduce the twelve-dimensional 
Dg+T^-equivariant bifurcation problems that enable us to determine the relative stability of simple 
hexagonal patterns, stripe patterns and certain rhombic and superlattice patterns. We make use of 
bifurcation results derived in |3^, which apply when there is a single critical wavenumber kc, 
to demonstrate how the magnitude of the cross-coupling terms are pivotal in determining pattern 
stability. As before, we consider a stroboscopic map, but now restrict analysis to patterns that are 
doubly-periodic on some hexagonal lattice. For instance, the free surface height takes the form 

h{-K,pT)= hMe"^"^'^'^'""^^^^"" + c.c. (11) 

meZ2 

at time t = pT, where ki,k2 €E generate a hexagonal dual lattice (|ki| = |k2| and ki • k2 = 
-i|kip); see Figure |[ 

The twelve-dimensional irreducible representations of Dg+T'^ apply to the bifurcation problem 
when there are twelve integer pairs (mi, 7712) in ( ITTI ) such that |miki + m2k2| — kc, where k^ is 
the critical wavenumber of the instability at the bifurcation point. See Figure ^ for an example. 
Following p3| we will associate with each twelve-dimensional irreducible representation an integer 
pair (711,712); in particular iii and 712 are co-prime, tii > 712 > 7ii/2 > 0, and tii + 712 is not a 
multiple of 3. The neutral modes that span the center eigenspace at the bifurcation point take the 
form 

{zi e'^i-" + Z2 e'^^-" + Z3 e'K^'-" + Z4 e*^^'" + Z5 e'^-^'" + e*^"" + c.c.|z, G C}, (12) 
where 

Ki = niki + 77.2k2, K4 = riiki + (rii - 7i2)k2, 

K2 = (-711 + 7l2)kl - 71lk2, K5 = -7l2kl - 71lk2, (13) 

K3 = -?i2ki + (ill — 7i2)k2, Ke = (712 — iii)ki + 7i2k2 . 

Note that ±Ki, ±K2, iKs point to the vertices of a hexagon, as do ±K4, ±K5, iKe, and that the 
two hexagons are rotated relative to each other by an angle 9h G (0, f ) indicated in Figure ^. This 
angle is related to (rii, 712) by 

Tl? + 2711112 - 2n2 

cos(6l,0 = — (14) 

2(71J - 711H2 + 712) 
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• • Kg • 

• • • • 

• • • 



Figure 3: Hexagonal k-space lattice, with critical circle of radius kc superimposed. In this example 
(ni, rt2) = (3, 2) in (p^, and the critical circle intersects twelve points that lie at the vertices of two 
hexagons rotated by 9h relative to each other. 



Also note that the ratio of lengthscales for superpatterns depends on (711,712). Specifically, |ki| 
determines the larger periodicity scale of the superpatterns, while |Kj| — kc determines the smaller 
lengthscale associated with the instability; thus the lengthscale ratio is 



|K,|/|ki| 



if — 7ii7l2 + 



(15) 



The example of Figure || corresponds to (711,772) — (3,2), for which 9h ~ 22° in ( |l^ and the 
lengthscale ratio (|lj) is the smallest associated with a hexagonal lattice, namely Vl. These are the 
angle and lengthscale ratio that apply to the experimental superlattice pattern reproduced from ||l[ 
in Figure 

The general form of the twelve-dimensional Dg+T^-equivariant mappings are derived in ]3^ . 
Through cubic order in Zj , they take the form 



Z2 
Z3 
Z4 

Z5 

z& - 



62122^ + 62^3^ + 6412:4'^ 

62|zi|^ + 62|23|^ + 64125 

62|zi|2 +62IZ2P +64IZ6 
62|25P +62I26P +64IZ1 



62|Z4|^ +62126^ +6412:2 



h\zb? + h\zii\^)zi 
h\z<i 

j23 



65IZ4 
65 1^3 
65IZ1 



'+b,\z,\' 



> (t{{1 + \)zi + ez2'z3 + (6i|zi 

> (t{{1 + \)z2 + fizizz + (6i|z2 

> a{{l + A)z3 + eziZ2 + (6i|z3 

> (7^(1 + A)z4 + ez5Z6 + (6i|z4 
y cr{{l + A)z5 + ez4Z6 + (6i|z5 

a[{l + A)z6 + ez4Z5 + (6i|z6p + 62|z4|^ + 62^5!^ + bi\zz? + 65|z2|^ + 66|zip)z6) , 



(16) 



2 +66|Z2P)Z4 
'+66|Z3|')Z5 



where A measures the distance from the critical excitation amplitude, and cr = +1(— 1) in the 
case of (sub)harmonic instability. All nonlinear coefficients are real. If ct — —1 then a normal 
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form transformation removes all even terms on the right-hand-side of (|T^) and hence e = 0. The 
dependence of the general equivariant bifurcation problem on (ni, 712) does not appear until higher 
than cubic order in its Taylor expansion [ p4| . 

We now recall some basic results pertaining to the bifurcation problem (p^. In the a — +1 case 
the equivariant branching lemma pl| ensures the existence of harmonic wave solution branches in 
the form of stripes, simple hexagons, rhombs, and super hexagons p3[ |. A primary solution branch 
with submaximal isotropy, named super triangles, was also shown to exist in [Q. See Figure |l|a 
for an example of this pattern. Table |l| gives the general form of these solutions, along with their 
branching and stability assignments. The general bifurcation results in the case that cr = — 1 can 
be found in |Q; this bifurcation problem differs from the harmonic case in that it possesses an 
additional Z2 normal form symmetry. The equivariant branching lemma then ensures existence of 
five additional solution branches to those listed in Table . 

The generic presence of a quadratic term in (|l^) for the harmonic case renders all of the solutions 
in Table 1^ unstable at bifurcation. Hence the transition from the flat state to the patterned harmonic 
wave state is expected to be hysteretic. In order to capture stable weakly nonlinear solutions, we 
must focus our analysis on the unfolding of the degenerate bifurcation problem e = 0. Note that 
when e = the stability of simple and super hexagons/triangles is not determined at cubic order 
since the phases (f>j of solutions zj = rje'"^^ to (^6|) are then arbitrary. Even in the case of < |e| <C 1 
the relative stability of super hexagons and super triangles depends on terms that are at least fifth 
order. However, we may use the cubic truncation to determine that one (and only one) of these 
two solutions is stable. The higher order terms are only needed to determine whether it is the 
hexagonal or triangular superpattern |Q. 

When < |e| ^ 1, it follows from Table |l| that a necessary condition for one of the superpatterns 
to be stable over some range of A values near onset is for 

&i + 2&2 < + + < . (17) 

The combination bi + 2^2 is independent of the lattice angle 9h in (pH); it is computed from 
a hydrodynamic model of the two-frequency Faraday problem in the Appendix by considering 
bifurcation to simple hexagons. In contrast, the combination 64 + 65 + &g depends on 9^ and is 
computed in the appendix from the hydrodynamic equations by considering the rhombic lattice 
bifurcation problem (^0|). Specifically, the cross-coupling coefficients 64,65,66 are 

64 = /?(0h), 65=/3(0„ + y), 66-/5(^?„-y), (18) 

where 9h is the angle between Ki and K4 given by ([T^). (The function l3{9) may be extended 
from 9 € (0, f ] to angles 9 € (0, 27r) using P{9) = P{-9) ^ (3{9 + tt), identities that follow from the 
symmetries of the rhombic lattice bifurcation problem.) 

The inequality ( p7| ) will be satisfied (if at all) only for those 9^ values where I64 + 65 + 65! is 
small compared to |6i + 262]. Moreover, if 61 — 62 < in addition to (p^), then simple hexagons 
become unstable on a given hexagonal lattice when 

€^(^4 + 65 + 66) 

(61 + 262 - 64 - 65 - 66)2 • ^ ^ 

If 64 + 65 + 65 < for all 9h, then simple hexagons first lose stability with increasing A to a 
perturbation in the direction of a superpattern for that value of 9h that minimizes [64 + 65 + 66|. 
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Table 1: Branching equations and stability assignments for the harmonic case (cr = +1); e, 61, . . . , 
are coefficients in the bifurcation equations (p^. A solution is stable if all quantities in the right 
column are nee ative. See [pj pj R for more details. 



Planform and branching equation 


old Dill ly 


Stripes 
z = (x, 0,0, 0,0,0) 

U — AX -f- OiX -\- LyyX j 


sgn(6i), 

sgn(ex + (&2 - bi)x'^), sgn(-ex + (62 - ^1)2:^), 
sgn(&4 — 61), sgn(65 — 61), sgn(66 — ^i) 


Simple Hexagons 

= Ax + ea;2 + {bi + 2b2)x^ + 0{x^) 


sgn(ex + 2(61 + 2b2)x^), sgn(-ex + (61 - b2)x^) 
sgni^— ex -|- ^04 + 05 + Og — Oi — Z02)x ) 
sgn(-ex + 0{x^)) 


Rhombs (Rh4) 
z = (x,0,0,x,0,0) 
= Xx + {bi+bi)x^ + 0{x^) 


sgn(6i+64), sgn(6i--64), sgn(Ci), sgn(C2), 
where Ci + C2 = (-26i - 264 + 2&2 + 65 + be)x^, 
C1C2 = -e^x^ + {bi +bi-b2- 65) (61 + &4 - &2 - bQ)x^ 


iXiioniDb 1^x1115 j 
z = (x,0,0,0,x,0) 


same as Rh4 with 64 65 


Rhombs (Rhg) 


same as Rh4 with 64 ^ foe 


Super Hexagons 

Z — ^ ijj ^ Xf ^ tXj ^ 'X ^ tXj ^ 'X ^ 

= Ax + ex^ + (61 + 2b2)x^ 
+ (&4 + ^5 + &6)a;^ + 0{x^) 


sgn(ex + 2(61 + 2&2 + &4 + 65 + fee)^^) 
sgn(ex + 2(61 + 262 - 64 - 65 - b(-)x^) 
sgn(-ex + e'(x3)), sgn(Ci), sgn(C2), 
where (1 + C2 = -4ex + 4(6i - 62)2;^, 
C1C2 = 4(ex - (61 - 62)x2)2 

-2((64 - 65)' + (64 - fee)" + (^5 - bQ)^))x^ 
sgn(C3), where Cs = 0(x2("i-i)) 


Super Triangles 

z = {z,z,z,z,z,z), 
z^xe''f',ij^O,Tr,... 


Same as super hexagons 
except Cs ^ -Ca 
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k k 

Figure 4: Neutral stability curves computed from ( ^o|) linearized about /i = $ = 0. Floquet 
multipliers of +1 (—1) are indicated by solid (dashed) lines, (a) m/n = 2/3, (p = 0° , x — Xc = 66.6°, 
To = 0.53, Go = 0.47 and 7 = 0.09 in (b) m/n = 6/7, = 0°, x = Xc = 53.0°, Tq = 7.5, 

Go = 1.5 and 7 = 0.08. 



If 64 + 65 + 65 > for any 9h, then small amplitude simple hexagons are unstable when A > 0. 
Thus we expect the stability properties of superpatterns and simple hexagons to be affected by the 
presence of a weakly damped harmonic mode when 9h or 6fi ± 27r/3 is near 9r (or t: ± 9r), the 
resonant triad angle, since it is in this situation that one of the cross-coupling coefficients 64, 65 or 
fog may suddenly change in magnitude. 



3 Results 

This section shows explicitly the role of resonant triads and weakly damped harmonic modes in 
the pattern selection problem for two-frequency forced Faraday waves. We examine how the cubic 
nonlinear coefficients in (|l6|), for the Zhang- Vinals hydrodynamic equations vary as a function of 
9h, the lattice angle and explain how this can be related to 9r, the resonant triad angle. The 
details of the computation of the coefhcients are relegated to the Appendix. We focus on two 
examples, involving forcing frequency ratios m/n = 2/3 and 6/7. The 2/3 case demonstrates the 
basic difference between the pattern selection problems for subharmonic and harmonic instabilities 
near the bicritical point. Our investigation also reveals a fundamental difference between harmonic 
wave pattern selection in the 2/3 and 6/7 cases, due to the presence of additional harmonic wave 
resonance tongues for the higher 6/7 forcing frequencies; see Figure ^. 

3.1 The Zhang-Vinals Hydrodynamic Equations 

The quadratic and cubic nonlinear coefhcients in the hexagonal bifurcation problem ( |l6| ) are com- 
puted in the appendix from a model of the two-frequency Faraday problem derived by Zhang and 
Vihals Q from the Navier-Stokes equations. Their equations, which apply to weakly damped, 
small amplitude surface waves on a semi-infinite layer of fluid, describe the evolution of the surface 
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height h{x,T) and surface velocity potential $(x, t). Specifically, 

drh = j\/^h + V^-\/ ■ {h\/^) + ^\/^{h^V<i>)-V{hV^)+V[hV{hV^) + ^h^\/^<^>] 

dr<p = 7v2$ + roV2/i-G(r)/t + i(i5$)2-i(v$)2-(i5$)[w2$ + p(/iP$)] (20) 
-iroV-((vM(v;i)2), 

where I? is a nonlocal operator that multiplies each Fourier component of a field by its wave 
number, i.e. Ve^^'^ = |k|e*'' ''. Here time has been scaled by uj so that the (non-dimensionalized) 
two-frequency acceleration is 

G(t) = Go — /(cos(x) cos(TOr) + sin(x) cos(nT + (f)j). (21) 

The damping number (7), capillarity number (Fq), gravity number (Go), and dimensionless accel- 
eration (/) are related to the forcing function (|^) and the fluid parameters by 

7= : l0 = J, Go = ^, / = 5-. [Z^) 

Here v is the kinematic viscosity, F is the surface tension, p is the fluid density, and the wave 
number k^ is chosen to satisfy the dispersion relation 

5ofco + — - . (23) 



p V 2 



3.2 Example 1: m/n=2/3 



This example demonstrates a result of the general normal form analysis of Section 2.2, namely 
that proximity to the subharmonic/harmonic bicritical point will strongly influence the pattern 
selection problem for subharmonic waves, but not for harmonic waves. Specifically, we examine 
the cross-coupling coefficient in ( |l0| ) as a function of the angle d for onset of both harmonic 
and subharmonic waves near the bicritical point. We show that only in the subharmonic case does 
1/3(6') I become large at the resonant angle 6',. in (||). 

As described in Section the primary instability changes from harmonic (Floquet multiplier 
+1) to subharmonic (Floquet multiplier —1) as % in ( |2l| ) is increased through the bicritical point Xc- 
This transition is determined from the linear hydrodynamic problem, which for the Zhang- Vifials 
model ( pO| ) takes the form of a damped Mathieu equation for each Fourier mode h = /ifc(r)e*'^^: 

h!l + llk^h!^. + (ci^k^ + ^X)hk = /cfc[cos(x)cos(mT) + sin(x) cos(nT)] /i^. (24) 

Here the natural frequency Sl^ satisfies the dispersion relation ri| ~ Ggk + Ygk^. A numerically- 
computed neutral curve f{k) for m/n = 2/3 forcing and X — Xc = 66.6° is given in Figure |4^. The 
other parameters of this example are = 0°, Fo = 0.53, Go = 0.47 and 7 = 0.09. 

We now vary x near Xc, holding all other parameters fixed, and examine the rhombic lattice 
cross-coupling coefficient P{9) in (p^ for onset subharmonic/harmonic waves, as appropriate. We 



13 



have scaled the ampUtudes vi and V2 in (|l^) so that a = — 1. We note that in the harmonic 
case (3 diverges as 9 ^ 60°, i.e. when the rhombic lattice approaches the hexagonal one and 
there is an additional mode associated with the center manifold dynamics. This is in contrast 
to the subharmonic case, for which there is a normal form symmetry that ensures existence of a 
dynamically invariant subspace spanned by a pair of subharmonic modes separated by 60°. Thus 
in the subharmonic case /3 remains finite at 6* = 60°. 

For X > Xc the primary instability is to subharmonic waves. For instance, for x = 66.7° the 
minimum of the neutral curve occurs at wavcnumber /cc.s = 1.415 with forcing amplitude fc = 0.842, 
and is associated with a Floquet multiplier cr = — 1. The nearly critical harmonic resonance tongue 
has its minimum at {k, /) — (0.962, 0.846). In this case, there is a spatio-temporally resonant triad 
comprised of the weakly damped harmonic mode and, from (0) , two subharmonic modes separated 



by 9r = 39.9°. It follows from our general analysis of Section 2.2 that P{9) will be large in magnitude 
for 9 near 9r. Figure ^ shows P{9) for this case, and indeed, the nonlinear coefficient exhibits a 
large dip centered at 6* = 6*,- — 39.9°. At this angle, \P{9)\ takes on its largest value. Similar 
observations have been made by Zhang and Viiials [23| for forcing frequencies in ratio m/n = 1/2. 

In contrast, when x < Xc so that the first instability to occur with increasing / is harmonic, 
we find that the weakly damped subharmonic mode leaves no signature in the plot (3{9). For 
instance, for x = 66.5° the primary instability is to harmonic waves at wavenumber kc.h = 0.963 
and forcing amplitude = 0.841. The subhharmonic resonance tongue has a minimum at (fc, /) = 
(1.415,0.843). While there is a spatially resonant triad involving two critical harmonic modes, 
which by (|^) are separated by 9r — 85.7°, the triad of modes is not spatio-temporally resonant. 
Figure ||b shows the cross-coupling coefficient (3{9) for this case (with the region near 60° removed). 
As anticipated, there is no signature of the weakly damped subharmonic mode in the plot. Similar 
observations have been made by Silbcr and Skeldon in the setting of one-dimensional surface 
wave patterns. 

3.3 Example 2: m/n=6/7 

This example demonstrates a fundamental difference between harmonic wave pattern selection for 
low forcing frequencies {e.g., 2uj/3uj) and for high forcing frequencies {e.g., Qlo/Ilh). This difference 
is due to the presence of multiple harmonic resonance tongues in the neutral curve associated 
with the higher forcing frequencies; see Figure ^. In particular, these resonance tongues suggest the 
possibility that weakly damped harmonic modes may influence the harmonic wave pattern selection 
problem. This is in contrast to the m/n = 2/3 example of the previous section, for which only 
subharmonic wave pattern competition was affected by weakly damped harmonic waves. In this 
section we also demonstrate that the weakly damped harmonic modes may stabilize harmonic wave 
superpatterns at a lattice angle 9^ ~ 9r, due to a near cancellation of the two terms that contribute 



to P{9r) given by (m as described in Section 2.2 



We focus on bifurcation to harmonic waves for x = 52.4°, which is close to the bicritical value 
Xc = 53.0°. The remaining parameters are = 0°, To = 7.5, Go = 1.5 and 7 = 0.08. We note that 
while the forcing frequency ratio m/n — 6/7 coincides with that used in the experiments of KudroUi, 
et al. [Q, the remaining parameters do not coincide with the experiment. One problem with using 
the experimental parameters in the Zhang- Viiials equations is that the primary instability then 
moves to a subharmonic resonance tongue at very small wavenumber, i.e., the first resonance 
tongue of Figure ^b. This is because the Zhang- Vihals model does not accurately capture the 
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Figure 5: Cross coupling coefficients (3{6) in ( ]lO| ) computed in the Appendix from ( pO| ) for the case 
m/n = 2/3 and (f> — 0° in (pT|). The fluid parameters used are given in the caption of Figure 

(a) X — 66.7° > Xc when the bifurcation is to subharmonic waves. Note the dip at 6* = 9r = 39.9°. 

(b) X = 66.5° < Xc, when the bifurcation is to harmonic waves. Because the (nearly) critical modes 
are not in temporal resonance, P{9) shows no special structure a.t 9 = 9r = 85.7°. We have removed 
from this plot the region near 9 = 60°, where (}{9) diverges. 



(b) 

, I , L 



damping at small k that is due to finite depth effects. 

In this example we find two prominent features in the plot of the cross-coupling coefficient P{9) 
in Figure ||a: a large dip at 9 = 67.6° and a small spike at 9 ^ 22.2°. We now discuss the origin of 
these two features. 

The large dip around 9 = 67.6° is not a consequence of two-frequency forcing. Specifically, 
the dip remains in f3{9) even for purely 6uj forcing (i.e. in the limit x ~^ 0); c/. plots of P{9) 
in Figures ^ and ^ which are obtained with x = 52.4° and x = 0°i respectively. Thus this 
feature may be understood in the context of single frequency forcing, and has in fact already been 
investigated by Zhang and Viiials in that setting. Specifically, if x = 0° then the forcing period 
is T' = ^ = ^ and the primary instability is to subharmonic waves with period 2T'. A plot 
of the corresponding neutral curve is given in Figure with the primary harmonic resonance 
tongue from Figure IJb superimposed on it. In this single-frequency setting the feature at 67.6° 
is understood as being due to the damped harmonic mode around k = 1.7 in Figure ||d. Perhaps 
more relevant to this discussion is our observation that this feature, which leads to a large value of 
1^4 + 65 + is destabilizing for superpatterns. To see this, we refer to the discussion surrounding 
equation ( |l7|) and to Figure which shows that 

0>6i + 262>54 + 55 + 66 = [3{9h)+l3(^9h + Y)+P{0h-Y)^ 

= f3{9h)+p{l-9h)+p[^+9h) for J + 0;, « 67.6°. (25) 

In contrast the spike at 6' = 22.2° in Figure ^a minimizes [64 + 65 + ^61 at 9h ^ 22.2°, as shown 
in Figure |b. As we show below, this feature can lead to a stabilization of superpatterns and a 
destabilization of the simple hexagons. First we provide strong evidence that the spike is due to a 
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resonance between the primary harmonic instabihty [a — 1) and a weakly damped harmonic mode 
with a real Floquet multipher /x that is close to 1 (see (^ and (||) of Section 2.2). In order to 



show this we must first compute the Floquet multipliers /i(fc) at the critical forcing amplitude fc 
to determine the wavenumbers k at which /i w 1 at the onset of instability. 

We determine the Floquet multipliers /x(fc) f = fc = 1.552 numerically from the linear 
problem (p^. These are presented in Figure |^. We find that the multipliers are well approximated 
away from the two primary resonance tongues by considering the unforced problem (/ = in 
equation |2j), for which 

^^^g27rA±^ A± = -7/c2 ± if7fc . (26) 

Figures 0a and ^ show the magnitude ^ and the phase of the Floquet multipliers /i = ^e*''' both 
as computed numerically from (|2j) (solid line) and approximated by (^) (dotted line). Figure |^c 
shows the real part of the Floquet multipliers, ^cosV', versus wavenumber fc. The "bubbles" in this 
plot correspond to wavenumbers at which the Floquet multipliers are real (as opposed to a complex 
conjugate pair). Weakly damped harmonic modes are associated with bubbles near a Floquet 
multiplier of +1. Numerically we find that there are small bubbles of real Floquet multipliers 
whenever the phase "0 is a multiple of tt; this is demonstrated in Figure In particular, we find 
a bubble at wavenumber k = 0.383, with associated real Floquet multiplier /j — 0.93. This mode is 
weakly damped and forms a resonant triad with primary harmonic modes separated by Or — 22.2° . 
(Here fcc,/i = 0.997 for the primary instability, which corresponds to fc„ in (||), with fc^ = 0.393 
determined by the weakly damped harmonic mode.) Here we have focused on the wavenumbers 
associated with real Floquet multipliers near = +1 since weakly damped modes with complex 
Floquet multipliers do not form a spatio-temporally resonant triad with the primary harmonic 
modes. 

We now present some hexagonal lattice bifurcation results for the specific parameters of this 
example, which are given in Figure The computation of the quadratic and cubic coefficients in 
the bifurcation problem (|l6|) is described in the Appendix. We scale the amplitudes zj in ([l^) so 
that bi = —1, in which case we find that e = 0.00014 and 62 = —2.73. Thus we expect results of 



Section 2.3, which focused on the unfolding of the degenerate bifurcation problem e = 0, to apply. 

We find that simple hexagons, super hexagons and super triangles all bifurcate transcritically 
with the subcritical branch turning around in a saddle-node bifurcation. The stripes and rhombs 
solutions arise in supercritical pitchfork bifurcations. These claims are true for all lattice angles 
dh since the cubic coefficients 61, . . . , 65 in (|l^) are always negative; see Figure Moreover, we 
find that simple hexagons are always stabilized in a saddle-node bifurcation and that they do not 
lose stability until after they reach the supercritical regime A > 0. In contrast, super hexagons and 
super triangles are always unstable at A = 0, since at that point the sign of the second eigenvalue in 
Table is determined by sgn(6i + 262 — 364 — 365 — 366), which is positive for all 6 (see Figure ||b). 
Thus, as A is increased through 0, we expect a jump to finite amplitude simple hexagons as the 
other primary branches of ( |l6| ) are unstable. 

We find that simple hexagons eventually lose stability as A increases since the following two 
expressions from Table |^ change sign to positive (at least for some 9^) 

sgn(-e2;+ (61 -62)2:^), sgn{-ex + {b^ + + be - bi - 2b2)x^). (27) 

as the amplitude x of simple hexagons grows with A. The first quantity changes from negative to 
positive at A « 3.2 x 10~®. The second quantity changes sign with increasing A only for those values 
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Figure 6: (a) Cross-coupling coefRcient P{9) in ( p^ ) computed from ( pO| ) for the case m/n = 6/7, 
= 0° and X — 52.4° < in (pT|), and for fluid parameters given in the caption of Figure^, 
(b) Plots of bi + 262 (dashed) and 64 + 65 + be (solid) versus 6h- We note that 9h only takes on 
the discrete values satisfying (|^. (c) Cross-coupling coefficient f3{9) for 6lj forcing only; we have 
used the same parameters as in (a) except that now x = 0°. (d) Neutral curve for single frequency 
forcing. Floquet multipliers of +1 (—1) are indicated by solid (dashed) lines, and are computed 
relative to the period T' — 2tt/6. The primary harmonic resonance tongue from the two-frequency 
case of Figure ^ is superimposed as a dotted Hne. 
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Figure 7: Floquet multipliers fi = ^e^^ computed from (|2^ ) for the parameters used in Figure and 
for critical forcing amplitude f = fc = 1.552. (a) Magnitude ^, and (b) phase "0 ^'s. wavenumber 
k. The solid lines in (a) and (b) are computed numerically, while the the dotted lines are obtained 
by considering the unforced problem / = 0; see equation |2^. (c) Numerically computed real part 
^cosV' of the Floquet multipliers. The "bubbles" correspond to real-valued Floquet multipliers. 
The boxed region, shown blown up in (d), reveals a tiny "bubble" around k = 0.383, with real 
Floquet multiplier /i = 0.93. 
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Figure 8: (a) A value at which the first (dashed) and second (solid) eigenvalue expressions in 




turn positive versus the lattice angle 9h- Note that simple hexagons (H) first lose stability to 
perturbations in the super hexagon/ triangle (SH/ST) direction at Oh « 22.2°. (b) Schematic 
bifurcation diagram for the (711,712) = (3,2) lattice of Figure |3[ which corresponds to 6h = 21.8°. 
Stable (unstable) solutions are indicated by a solid (dotted) line. We do not show secondary 
branches or primary branches that are never stable. The stable rhombs solution (Rh) corresponds 
to one with an angle of 81.8°, which is the rhombs solution closest to 90° for this hexagonal lattice. 
The other two rhombs solutions are unstable. 



of 9h where 64 + ^5 + 65 ~ &i — 262 > 0, a condition which is met for 6^ > 11.5°. Figure ^ shows 
the value of A where the expressions of ( ^7|) change sign as a function of Oh. It follows that simple 
hexagons lose stability first on the lattice with angle 9h ~ 22.2°. This instability has an associated 
eigenvector in the direction of super hexagon/triangles, and at this value of A, super hexagons (or 
triangles) are stable. These results are summarized in Figure which shows part of the bifurcation 
diagram computed for the hexagonal lattice with (ni,n2) = (3,2), which corresponds to an angle 
dh — 21.8°. Note that when simple hexagons lose stability, both rhombs (Rh) and a superpattern 
are stable. Because the instability that first destabilizes the simple hexagons is in the direction of 
a superpattern with 9^ ~ 22.2°, we expect that the transition would be a hysteretic one involving 
the simple hexagons and a superpattern, at least in the absence of noise and other imperfections. 
We cannot determine whether the superpattern is hexagonal or triangular from our calculations, 
since this requires knowledge of fifth order terms in the bifurcation problem ]3^ . 



4 Conclusions 

In this paper we have examined the effect of spatio-temporally resonant triads on two-dimensional 
pattern selection in parametrically excited systems. Using a normal form transformation to enforce 
temporal symmetry and center manifold reduction, we have argued that weakly damped harmonic 
modes can strongly influence pattern selection by causing certain cubic cross-coupling coefficients 
in a twelve-dimensional Dg-|-T^-equivariant bifurcation problem to suddenly vary in magnitude for 
certain lattice angles 9h ■ This suggests an important consideration in choosing one over another of 
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the countable set of twelve-dimensional representations relevant to hexagonal bifurcation problems. 
Weakly damped subharmonic modes, on the other hand, do not have such an effect. 

Our general analysis applies to any parametrically excited pattern forming system, but in par- 
ticular is relevant to the interpretation of many recent experiments on two-frequency forced Faraday 
waves. In such experiments, a bicritical point exists where subharmonic and harmonic instabilities 
are simultaneously excited. On one side of the bicritical point, a subharmonic mode is excited and 
there is a weakly damped harmonic mode, while on the other, it is the harmonic mode which is 
excited and the subharmonic mode which is weakly damped. We showed that this weakly damped 
subharmonic mode does not influence the harmonic wave pattern selection problem. 

We have derived the quadratic and cubic coefficients in the rhombic and hexagonal bifurcation 
equations describing the onset of patterns from the hydrodynamic equations of Zhang and Vihals. 
We presented results for two different sets of parameters. In the first case, the two forcing frequencies 
are in the ratio 2/3 and the modes near the bicritical point are the only ones of relevance. As 
expected from our normal form analysis, for subharmonic waves, the weakly damped harmonic mode 
affects the cross-coupling coefficients, while for harmonic waves, the weakly damped subharmonic 
mode had no effect. 

In the second case of 6/7 forcing we have shown that, in addition to the modes near the bicritical 
point, there are other harmonic modes that are important. These modes are not close to onset in 
the sense that they only become critical at a much higher value of the excitation amplitude, but 
are weakly damped and thus must be taken into account. We demonstrate that they can have a 
stabilizing effect on superlattice patterns at a lattice angle approximately equal to the angle of the 
harmonic-harmonic resonance. This can occur if the contribution of these weakly damped modes 
to the nonlinear cross-coupling coefficient nearly cancels the other contributions to this term, and 
hence is a subtle effect that depends on certain details of the nonlinear problem, as well as the 
results of the linear analysis which identifies the near critical modes. For the parameters we have 
chosen, the onset pattern is simple hexagons, but upon a further increase of the forcing, there is an 
instability to a superlattice pattern associated with a hexagonal lattice with (rii, — (3, 2). 

The experiments of KudroUi, Pier and GoUub [Q found a superlattice pattern near the bicritical 
point which sits on a lattice with (ni,n2) = (3,2). The work in this paper suggests that the 
observation of this pattern could be explained by the interaction of the primary harmonic instability 
and weakly damped harmonic modes. However, the Zhang- Vifials equations are not valid in the 
parameter regime where this experiment was performed, and thus a study of the full hydrodynamic 
problem is necessary to confirm this conjecture. A complete study should also involve a more 
complete analysis of the codimension-2 bifurcation point and the associated dynamics, in the spirit 
of J.D. Crawford's early work on competing instabilities in the Faraday problem P7|. This would 



be of interest in light of recent two-frequency experiments by Arbell and Fineberg that show 
a variety of dynamic states near the bicritical point, which involve both critical modes. 

A Perturbation Theory 

Here we outline the computation of the coefficients in ([l^) and ( [l^ ) from the equations of Zhang and 
Vinals (pO|). A multiple-scale perturbation method is used to derive expressions for the coefficients 
which are then evaluated numerically using a pseudospectral approach. This follows closely the 
method described in p6[ for the onset of one-dimensional patterns and we refer the reader there 
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for further details. 

The coefficients can be derived by considering two different calculations, namely the bifurcation 
problem ( ]l6| ) restricted in turn to a rhombic and a simple hexagons subspace. 

A.l Rhombic lattice computation 

In order to compute the coefficient a and the cross-coupling coefficient f3{9) in ( p^ ) we seek solutions 
which are periodic on a rhombic lattice associated with an angle 9. We are thereby able to compute 
the coefficients bi, 64, 65, and feg in the bifurcation equations ( p^ ) since bi = a, b^ = P{9h), 
&5 = PiOh + 27r/3), and be = Pi9h - 27r/3). 

First we introduce a small parameter rj, such that 

H^,y,T) = ■qhi{x,y,T,T) + jfh2{x,y,T,T) (28) 

+ 7fh:i{x,y,T, T) + --- 
Hx,y,T) = 77$i (x, y, T, T) +772$2 (a;, y,r,T) 

+ r;3$3(2;,2/,r, T) + • • • , 

in ( pO| ) where 

T = v\ .f = fc + v\f2. (29) 

Here fc is the critical excitation amplitude. The terms in the expansion for h and $ may be written 
in the following separable Floquct-Fourier form: 

hi = [u.i(r)e''^'^^ + i«4(T)e*'=<=('=^+''2') +c.c.]pi(r) (30) 

$1 = [u;i(r)e''='=^ + u;4(T)e*'=<=('=^+''2') +c.c.]gi(T) 

h2 = K(r)e2^'=<=^ + w^(T)e2''='^(^^+''^)]p2,i(r) 
+ w;i(r)W4(T)e''=^«i-^)"-^^)p2,2(r) 
+ u;i(T)?«4(T)e*'==«i+^)"+^^V2,3(T) + c.c. 

$2 = K(T)e2*^-=" + u;2(r)e2*'==(""+"^)](724(T) 

+ u;i(T)?«4(T)e*'=-«i+=)"+^^)g2,3(T) + c.c. 

where c = cos0, s = sin 9, and 9 is not a multiple of ^. Here pi and gi are real 27r-periodic 
functions of the fast time r in the case of harmonic waves; in the case of subharmonic waves they 
are 47r-periodic in r. Additionally, p2.r and (72, r {r ~ 1, 2, 3) are real 27r-periodic functions of r. 
The wave number kc is associated with the onset unstable mode. 

At 0{r]) we recover the linear problem which determines kc and fc, as well as the functions pi, 
qi to within a multiplicative constant. At O^q^), equations are found which allow us to solve for 
the functions p2,r and q2,r- Finally, at 0{rf'), we apply a solvability condition to ensure that a 
periodic solution exists. This condition leads to the amplitude equations 

af2Wi+ A\wi\^wi+B{9)\wi\^wi (31) 

cthwi + A|w4pu;4 + B{9)\wi_\^wa, 



,dw4 

'It 
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where 



1 

^ J (P'l +lklpi)pi dr 



a 



A = 



27r 

k 

U2 



(32) 



[cos(x) cos(mT) + sin(x) cos(nr + dr 



^ y [-fcc(p?gi)' - iklplqi - 2(gip2,i)' " 27fcc'?iP2,i 



B{e) 



in 



r r 

/ (1 - c - \/2 - 2c)[(pig2,2)' + 7fccPi92,2 - kcqiq2,2\ 
JO '- 

+ (1 + c - \/2 + 2c)[(pig2,3)' + 7fccPi92,3 - fccgig2,3] 

- (1 - C)[(p2,2gi)' + 7^cP2.2'?l] - (1 + C)[(p2,3«l)' + 7fccP2,3'7l] 

- (6 - 2^/2~2^- 2V2T2^)[fcc(p?gi)' + 7fccP?'Zi " klPiQi] 
+ Toi3c'^ + s'^)klpl)]pi dr. 



In the above, a prime denotes differentiation with respect to r and pi is the equivalent of pi for 
the adjoint problem at 0{ri). The amplitude equations (|3l| ) may be re-scaled and then comparison 
with the map (M) yields 



a — bi — sgn(AQ!), f3{9) = sgn(AQ!) 



m 

A 



(33) 



A. 2 Hexagonal lattice computation 



Similarly, we compute the coefficients e and 62 in the bifurcation equations by seeking solutions 
in the form of simple hexagons. Here we use a three-timing perturbation method, writing the 
solution as 



where 
and 



hi 
h2 



Hx,y,T) = T]hi{x,y,T,Ti,T2) + ■n'^h2{x,y,T,Ti,T2) 

+ T]^h3{x,y,T,Ti,T2) H 

<f{x,y,T) = r]<S>i{x,y,T,Ti,T2) + v^^2{x,y,T,Ti,T2) 

+ r;3$3(x,y,T,Ti,T2) + ---, 

Ti = ryr, T2 = rfr, 

Wi(ri,T2)pi(r)[e*'=-^ -I- e*''-<=(-5=^+'?2') -t- e*'=^<=(-4^-'??^) + c.c.; 
i«i(ri,T2)gi(r)[e^'=-^ -f e''^=(-^^+'5^J^) + e^'^-'f-^^-^^) + c.c] 
wl{Ti,T2){p2AT)[e'^-'^'' + e'^-'^-'-'+^y^ + e'^-^-'^'^y^ + c.c 



(34) 



(35) 
(36) 
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- c.c] 




_|_ ^ikaVSy 


gikai^x- 




C.C.]} 


wl{T,,T2){q2ArW''^''' + e"'^^- 


x+VSy) _j_ 








ix+^y) ^ 






- C.c] 


+ g2,3(r)[e*'=^(^"^'?^) 




_|_ g^fccd^- 




C.C.]} 



As with the rhombic case, pi, qi, p2,r and 52, r are real. Additionally, we take the amplitude 
wi(Ti,T2) to be real. 

For the harmonic case, at 0{rf ) the solvability condition, 



must be satisfied, where 5 is given by The quadratic coefficient is 



f^o = ^ J [-{PiqiY - iklpiqi + -kcql]pi dr 



(37) 



(38) 



There is no solvability condition for subharmonic waves at 0{t]^), reflecting the fact that there are 
no even terms in the amplitude equations ( p^ for this case. 

At Oiji^), we again apply a solvability condition to ensure that a periodic solution exists. This 
conditions leads to the amplitude equation 

.dwi 



a ±2 



The coefficients 5, a, and A are given by (p2), and 



47Z ^ 2 

(2 ~ ^/3)fcc [(Pl(72,3)' + 7^cPl'72,3 - fccgi<?2,3] 

(2V3 - 4)A;3[(p2gi)' + „ k^p^ql] 

3 

T^fcc [(P2.3'7l)' + 7fccP2,39l - Tofc^p?] 



(39) 



(40) 



:^c[(Pi'?2,2)' + lklpiq2a + (^2,2(71)' + lklp2.,2li - fcc<7ig2,: 



/3o r, 2 , /^O , „ , I r, ; 2 1 

— [fc^Pigi + -jPi + 2p2,2 + 27«:c-P2,2j 



Pi dr. 



By rescaling r]Wj{Ti, T2) Wj{T) and rf'af2 a/2, we obtain the reconstituted hexagonal bifur- 
cation equation 



= a/291 + Pawl + (A + 2(32)wl 
Finally, after rescaling as for the rhombic case, and comparing (^l]) to ( |l6| ) we find that 

e = sgn a) , 62 = sgn(Aa) — . 



(41) 



(42) 
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